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About these notes 
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Preamble 

Classical physics - Newton's laws of mechanics and the theory of electromagnetism developed 
in the 19th century - fails at the atomic scale. From the beginning of the 20th century, mounting 
experimental evidence pointed towards the existence of a radically different theory of physics 
that governs the properties of atoms and their constituent particles. The deterministic and 
continuous nature of classical physics was replaced by a new quantum theory, with probability 
and discreteness at its heart. The theory of quantum mechanics developed in the 20th century 
not only successfully describes the structure of atoms and molecules, but also nuclear physics, 
particle physics (such as in particle accelerators like the LHC), chemistry (such as chemical 
bonding), the structure of solids, superconductors, etc. Quantum theory also underpinned many 
important technological advances in the 20th century, such as the laser, the microchip (hence 
computers and mobile phones), and the electron microscope. Future applications may include 
quantum cryptography and the quantum computer. 
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In this course we begin with an overview of some key physical ideas and formulae. These 
developed from experiments that demonstrate the failure of classical physics at the atomic scale. 
An important concept here is wave-particle duality. This, together with some intuition from 
classical physics, will lead us to the Schrddinger equation that governs such matter-waves. We 
discuss general properties of the Schrodinger equation, and interpret the wave function in terms 
of a probability distribution. We also study in detail some of the simplest, and most important, 
solutions: a particle confined to a box, the harmonic oscillator (which universally describes small 
oscillations of any quantum system), and finally the hydrogen atom. In particular, we will derive 
the observed emission/absorption spectrum of the hydrogen atom, a computation that was in 
Schrodinger's original 1926 paper. 

Classical particles and waves 

Before starting the course proper, we begin with a brief review of classical point particles and 
waves. This is to remind you of some concepts and formulae learned at school or in Mods, 
and to set notation used later in the text. A summary of vectors, and in particular differential 
operators such as V = grad and the Laplacian V^, is contained in appendix A. 

0.1 Point particles 

A point particle is an idealized object that, at a given instant of time, is located at a point in 
space. Considering a time interval [to,ti] C M and modelling space by M^, the particle's motion 
is described by r : [to,ti] — >■ R^. This is governed by Newton's second law: if the particle has 
constant mass m and is acted on by a force F, then r{t) obeys 

d^r 

F = m^. (0.1) 

In principle here F = F(r, t), although in this course we shall only consider conservative forces 
where F = —W, for some function V = V{r) called the potential. In particular, notice that V 
is independent of time (static). In this case, the total energy of the particle 

E = ^mlfp + F (0.2) 

is conserved, where we shall sometimes denote r = dr/dt. That is, E is independent of time: 
dE dr d^r dr dr / d^r 



, - m- — + VF- — = — • m— ^-F =0. (0.3) 

dt dt dt^ dt dt \ dt^ J ^ ^ 

Here the particle has kinetic energy ^m|rp, potential energy V = V{r(t)), and momentum 
p = mr. We may also write the kinetic energy as 

|p|2 

£^kinetic = ^ . (0.4) 

A particle that is subjected to no forces is called a free particle. 
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* The above formulae are true (to a very good approximation) for a non-relativistic point 
particle, which means its speed u = |r| is much less than the speed of light, u <^ c = speed 
of light. For example, a massless particle, with m = 0, necessarily moves at speed c and 
has energy £^ = c|p|, where p is its momentum (see relativity course next year). 

0.2 Waves 

Recall the classical wave equation 
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= V^^, (0.5) 



where v is the constant speed of the wave. This linear equation governs, for example, the 
propagation of sound or light. As basic solution we have the complex plane wave 

(j){r,t) = ^ exp [i(k • r - wi)] , (0.6) 

where we have the constant wave vector k, angular frequency u, and (complex) amplitude A. 
Substituting into (0.5) gives —cj^/v'^ = — |kp, or equivalently the relation 

" = M- ^'-'^ 

Both the real and imaginary parts of 4> separately satisfy the wave equation, giving real solutions 
that arc linear combinations of sines and cosines in (k • r — ojt). In fact it is a result of Fourier 
analysis that every solution to the wave equation (0.5) is a linear combination (in general 
involving an integral) of these plane waves. 

The wave frequency isv = uj/2'!r, while the wavelength is A = 27r/|k|, so that (0.7) is equivalent 
to saying that the wave speed is v = vX = frequency x wavelength. 

1 Physical background and wave-particle duality 
1.1 The photoelectric effect and Einstein-Planck relation 

In the mid 19th century Maxwell successfully described light as a wave propagating in the 
electromagnetic field. We shall not need to know anything about electromagnetic theory in this 
course. Nevertheless, we note in passing that, in Maxwell's theory, light propagating through 
a vacuum is described by two vector fields E(r, t), B(r, t), called the electric and magnetic 
fields, respectively, each Cartesian component of which satisfies the wave equation (0.5), with 
V = c being the speed of light in vacuum. This theory of clcctromagnetism unified the theories 
of electricity, magnetism and radiation, and explained wave-like properties of light such as 
reflection, polarization and diffraction. However, by the beginning of the 20th century it was 
becoming clear that Maxwell's theory could not explain experiments at the atomic scale. 

A clear and simple experiment that demonstrates this is the photoelectric effect, shown in 
Figure 1. Light of angular frequency uj is shone on a metal plate in vacuum. Electrons e~, 
which are only weakly bound to the metal, are emitted from the surface. One measures their 
kinetic energy and discovers the formula 

£^kinetic = -^0 + ^ • (1-1) 
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Figure 1: To observe the photoelectric effect, light is shone on a metal plate in vacuum. Electrons 
e~ are emitted from the surface of the metal, and their kinetic energy .Ekinetic is measured. 

• £^0 > is a constant energy which depends only on the particular metal used. 

• The constant of proportionality h (usually read as "h bar" ) is a constant of Nature that is 
fundamental to Quantum Theory. Prom (1.1) we see that it has dimensions [h] = energy 
X time, or equivalently [H] = M{LT~^)'^ x T = MLT~^ x L = dimensions of angular 
momentum, where M, L and T denote dimensions of mass, length and time, respectively. 
Numerically, ^ 2± 1.05 x 10~^^ Js. The combination 2Trh = h is called Planck's constant, 
while H is sometimes referred to as the reduced Planck's constant. 

• If the angular frequency u < Eq/H, no e~ are emitted. 

• The formula (1.1) is independent of the intensity (brightness) of the light, but as the latter 
is increased the number of e~ emitted increases. 

The classical theory of Maxwell does not explain these observations. Instead Einstein made 
the following remarkable hypothesis: 

Light of angular frequency to exists in small packets, or "quanta" , of energy 

E = hw (Einstein-Planck relation) . (1-2) 

These packets of light are known as photons, and are massless particles (m = 0) that travel at 
the speed of light c (c/. the starred remark at the end of section 0.1). In Einstein's interpretation 
of the photoelectric effect, each electron e" absorbs one photon of energy Tiu) (very occasionally 
more than one). Part of this photon energy goes into overcoming the "binding energy" Eq of 
the electron to the metal; the remainder is then converted into the observed kinetic energy of e~ 
when it is emitted. Increasing the intensity of the light simply increases the number of photons. 
Equation (1.2) is called the Einstein-Planck relation.^ 

^* Planck had introduced the constant h = 2'kH in earher work on radiation emitted from certain idealized 
hot objects, called black bodies. 
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1.2 Emission/absorption spectra of atoms 

Atoms emit and absorb light at very particular frequencies. The simplest atom is the hydro- 
gen atom, which we shall study in more detail at the end of this course. For hydrogen these 
frequencies were discovered experimentally in the 19th century, and are given by the formula 



Here ni < n2 are positive integers, and i2o — 1-10 x 10 m is Rydberg's constant, named after 
the discoverer of the empirical formula (1.3). 

Prom Einstein's description of light in terms of photons, this implies that a hydrogen atom 
emits and absorbs photons of particular energies hjjJni,n2- By conservation of energy, the en- 
ergy of the hydrogen atom itself must then be changing by these amounts when a photon is 
emitted/absorbed. This strongly suggests that the energies of the hydrogen atom must be given 



where n is an integer, so that e.g. for n\ < n2 the energy of the atom can change from 
to by emitting a photon of frequency L0ni,n2- That the hydrogen atom energies are indeed 
quantized in this way - that is, taking particular discrete values, rather than being continuous 
- has been confirmed in many other experiments. At the end of this course we will derive (1.4) 
theoretically. 

1.3 The double slit experiment 

The photoelectric effect implies that light, described classically by waves satisfying the wave 
equation (0.5) with speed v = c, is sometimes better described as a beam of particles, namely 
photons. 

Similarly, particles, such as the electron e~, can display wave-like characteristics. Perhaps 
the best experiment that demonstrates this is the famous double slit experiment, performed with 
electrons. This is shown in Figure 2. A beam of electrons is fired at a double slit configuration, 
with a detector screen on the other side. An electron hitting the screen appears as a bright spot, 
and over time one can plot this as a distribution. The latter exhibits a familiar diffraction pat- 
tern, similar to that seen in the corresponding experiment with the beam of electrons replaced 
by a beam of light. Such diffraction patterns are explained by the interference of waves: two 
waves that travel through each of the slits and arrive at the same point on the detector screen 
have travelled different distances. These waves then either constructively or destructively inter- 
fere with each other, depending on whether the difference in these distances is an even or odd 
number of wavelengths, respectively. This is perhaps familiar to those who have done A-level 
physics. 

A remarkable point here is that the diffraction pattern is still observed even when only a 
single electron is passing through the slits at a time. In fact this is the case in the Hitachi 




(1.3) 
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(1.4) 
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(a) Experimental set-up 
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(b) Actual electron hits in an experiment by a group 
at Hitachi (© Hitachi, Ltd): (i) 8 e", (ii) 270 e", 
(iii) 2,000 e", (iv) 160,000 e". 



Figure 2: The double slit experiment, performed with a beam of electrons. 



experimental results shown in Figure 2(b) (the time lapse up to picture (iv) is 20 minutes). This 
implies that the electrons are not interfering with each other to cause the diffraction pattern, 
but rather a single electron is behaving like a wave. Or, more precisely, the detected distribution 
of electron particles is characteristic of a wave passing between the slits. Notice that we may 
also interpret this distribution as a probability distribution for where any single electron will hit 
the screen. In this viewpoint, it is the probability that displays wave-like characteristics, while 
the electrons themselves are always detected on the screen as localized particles. These remarks 
are absolutely central to wave-particle duality. 

The double slit experiment (and variants of it) is extremely interesting and subtle, and we refer 
the interested reader to the references (especially the Feynman lectures) for further discussion 
of its role in understanding wave-particle duality. 

1.4 De Broglie's matter- waves 

The experiment we have just described suggests that particles, such as electrons, are also asso- 
ciated with waves. De Broglie made this more precise: 

A free particle of energy E and momentum p is associated with a wave of angular 
frequency co and wave vector k via 



E = huj , 

p = hk (de Broglie relations) . (1.5) 

Since the wavelength is A = 27r/|k|, we may also write the latter relation as A = 27r/i/|p| = h/\p\. 
De Broglie's insight was that these relations should apply to all particles, not just massless 
photons. In this context, E = hoj is usually referred to as a de Broglie relation, rather than the 
Einstein-Planck relation. 
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* We note that for a photon, the second relation in (1.5) is imphed by the first relation 
E = fujj. This follows from the starred comment at the end of section 0.1: for a photon 
E = c|p|, so that E = tujj together with c = w/|k| implies that |p| = fi|k|. In Special 
Relativity, which you can study next year, you'll learn that (£', p) and (a;, k) are both 
4-vectors, and indeed this was part of de Broglie's reasoning. 

2 Wave mechanics 

2.1 The Schrodinger equation 

De Broglie had hypothesized that particles, such as the electron e~, are associated with waves. 
Schrodinger set out to discover the equation that governs these matter-waves. He began by 
considering the plane wave, reviewed in section 0.2, 

*(r,t) = ^ exp [i(k • r - wt)] . (2.1) 

This of course satisfies the wave equation (0.5), with <j) replaced by the change of notation is 
meant to emphasize that we now wish to reinterpret this plane wave as a de Broglie matter-wave. 
Making use of the de Broglie relations (1.5), we notice 

ih— = huj'^ = E^ , 
dt 

-ihV^ = nk^- = p^r . (2.2) 

In particular, taking the divergence of the second equation implies that — ^^V^^* = |pp^', where 
is the Laplacian. 

The de Broglie relations apply to a free particle, for which the particle's energy is equal to 
its kinetic energy (since F = the potential is constant, and we take this constant to be zero). 
If the particle has mass m, then (0.4) relates E = ^^kinetic = |Pp/2?Ti. Using the de Broglie 
relations (1.5) then implies 

CO = - = . (2.3) 

Compare this to the relation (0.7). Putting everything together, we have 

ih^ = E^ = ^-^^ = -^V^* . (2.4) 
dt 2m 2m ^ ^ 

Thus we may associate to a free particle of mass m a plane wave ^'(r,t), using the de Broglie 

relations, which then satisfies the equation 

= -<!^V^*. (2.5) 
dt 2m ^ ^ 

This is essentially expressing the relation (2.3), which is simply the relation between energy and 

momentum for a free particle. 

More generally, a particle of mass m moving in a potential V = V{r) has energy 

is = UJ + y. (2.6) 

This led Schrodinger to 
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Definition / postulate A single, non-relativistic particle of mass m moving in a potential 
V{r) is described by a wave function '^{r,t) that is governed by the Schrddinger equation 



The wave function ^{r,t) is precisely de Broglie's matter-wave. 

Let us make some remarks about what wc have done above, which is deceptively straightfor- 
ward, and about the Schrodinger equation itself: 

1. It is important to realize that we have not derived the Schrodinger equation, in any rigorous 
sense. In particular, in the last step we have taken formulae that apply to free particles, 
that have been suggested by a combination of experiments and theoretical arguments, and 
extrapolated this, using the classical formula (2.6) for energy, to an equation governing 
the matter-wave of a particle moving in a general potential. It turns out this equation is 
indeed correct, but nothing in the rather naive argument we gave really guarantees this. 
The real test of the Schrodinger equation is that it agrees with experiments. 

2. The Schrodinger equation is a linear partial differential equation for a complex-valued 
function ^'(r,t). Thus if *i,^'2 are solutions, then so is q;i*i -|- a2^2, for any complex 
constants ai,a2 € C. This implies that solutions form a (usually infinite-dimensional) 
vector space over C. It is precisely this superposition of wave functions that leads to 
interference effects, as in the double slit experiment with electrons. Notice that, in contrast 
with the classical wave equation (0.5), the Schrodinger equation (2.7) is complex, due to 
the i = on the left hand side. 



3. Although we began our exposition by discussing photons, it is important to remark that the 
photon is a massless, relativistic particle, and as such is not governed by the Schrodinger 
equation. The quantum theory of photons is a much more involved theory, known as 
quantum electrodynamics, that requires a thorough understanding of both classical elec- 
tromagnetism and Special Relativity, as well as quantum ideas. We shall only refer to the 
photon again in the context of emission/absorption in atoms, for which we need only the 
Einstein-Planck relation. 

Before continuing to discuss some basic mathematical properties of the Schrodinger equation, 
and looking at our first example, let us pause to comment on the change of viewpoint that is 
already implicit in what we have said so far. Consider the classical problem of a point particle 
of mass m moving in a potential V. The dynamics is governed by Newton's law (0.1), with 
F = — W, the solutions of which give the particle's trajectory r{t). For given initial conditions, 
say the particle's position r and momentum p = mr at time t = to, one solves for the trajectory 
r{t), which gives the particle's location and momentum at any subsequent time. 

The corresponding quantum mechanical problem is very different. Given a quantum point 
particle of mass m moving in a potential V, we should instead solve the Schrodinger equation 
(2.7). We shall discuss the boundary conditions involved later, but notice immediately that the 



~dt 



2m 



(2.7) 
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result will be some complex- valued function *(r, t). You might immediately wonder what this 
function has to do with the particle's position at some time t. Again, we shall address this 
shortly. 



2.2 Stationciry states 

It is natural to seek separable solutions to the Schrodinger equation. Thus we write ^(r, 
V'(r)r(t), so that the Schrodinger equation (2.7) becomes 



(2.8) 



Since the left hand side depends only on t, while the right hand side depends only on r, both 
sides must be constant. If we call this constant E (anticipating that this will be the energy of 
the particle), then in particular we have 

i^^ = ET, (2.9) 

which immediately integrates to 

T{t) = e"'^*/'^ . (2.10) 

Here we have absorbed the overall multiplicative integration constant into ■0. The full wave 
function is thus 

*(r,t) = il;{r)e-'^^/^ . (2.11) 

That E is then indeed the energy of the particle follows from the de Broglie relation between 
energy and angular frequency for matter- waves: for this wave function the angular frequency is 
Lo = E/h, or equivalently E = huj. The function ij) then satisfies 

Definition The time-independent, or stationary state, Schrodinger equation for a particle of 
mass m and energy E moving in a potential V = V{r) is 

-—V'^ip + V^ = EtP . (2.12) 
2m 

The wave function ^'(r, i) = ^/;(r) e~'^*/'* is then called a stationary state wave function of 
energy E, although in a common abuse of language the function ijj is also often referred to as 
the stationary state wave function. 

2.3 One-dimensional equations 

Although ultimately we are interested in studying particles moving in the three spatial dimen- 
sions that we observe, it is often technically more straightforward to study the one- dimensional 
Schrodinger equation 



with corresponding stationary state equation 



2m dx^ 



and 



(2.14) 



(2.15) 



Here wc have replaced the Laplacian by the corresponding one-dimensional operator, which 
is simply 'd^ jdx? (see appendix A). Equations (2.13), (2.14) govern a particle propagating on 
the X-axis with potential V = V{x). Similar remarks apply in two dimensions. Although the 
one-dimensional equation looks somewhat unphysical, in fact sometimes a three-dimensional 
problem effectively reduces to a lower-dimensional Schrodinger equation; for example, due to 
symmetry reduction (see section 5), or because the particle is constrained to lie in some subspace. 

2.4 Pcirticle in a box 

Consider a particle inside a "box" on the x-axis. This means that the particle moves freely 
inside some interval [0, a] C M, but cannot leave this region. One can model this by a potential 
function V = V{x) that is zero inside the interval/box, and infinite outside: 



y(x) 



, < X < a , 

-|-oo , otherwise . 



(2.16) 



This is also sometimes referred to as the infinite square well potential. See Figure 3. 




V(x) = 




Figure 3: A particle in box [0, a] on the x-axis. 

Before discussing the quantum problem, let us briefly comment on the classical problem. Since 
V = inside the box, the classical particle moves at some constant velocity, or equivalently 
constant momentum p. Since the energy E = p'^/2m is conserved when the particle hits the 
edge of the box, after the collision p is replaced by —p and the particle heads towards the other 
edge of the box. Classically, notice that E may take any non-negative value. 

Now let us consider the quantum particle. The stationary state Schrodinger equation (2.14) 
inside the box is 

dV "^rnE 



dx^ 



(2.17) 
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for X G (0, a). We shall discuss boundary conditions more systematically in section 3.2, but here 
we note that since F = oo outside the box, the Schrodinger equation will make sense only if 
ijj = there. ^ When we come to discuss the physical meaning of the wave function, we shall see 
that this assertion is very well justified physically. If we also assume that ip is continuous, then 
we must solve (2.17) subject to the boundary condition V'(O) = V'('*) = 0- 
The general solution to (2.17) is 



' Acos^^x + Bsm^^x , 
A + Bx , 

A cosh ^^f^x + B sinh ^^f^x 



E>0 
E = 
E <0 



(2.18) 



In all cases the boundary condition ip{0) = implies ^ = 0. When E < 0, ip^a) = implies 
also B = 0. Thus the only solution is ^p = 0, which is always a physically meaningless solution 
to the Schrodinger equation (again, we shall discuss this more later). On the other hand, for 
E > the boundary condition 'tp{a) =0 implies that (either B = and V' = or) 



V2mE 



mr 
a 



for some integer n G Z. Thus the solutions 




< X < a , 
otherwise , 



(2.19) 



(2.20) 



are labelled by n. Notice that, without loss of generality, we may take n > 0. The first three 
wave functions, for ri = 1, 2, 3, are shown in Figure 4. 



B 






Figure 4: Wave functions for the first three states of a particle in a box. 'tpi{x) is the ground 
state wave function. 



The associated energies are, from (2.19), 



E = En = 



(2.21) 



2ma'^ 

We see immediately that the energy is quantized, i.e. it takes values in a discrete set, here 
labelled by a positive integer. This is in stark contrast with the energy of the classical particle. 



This is not a very rigorous statement. We can make the discussion rigorous by simply declaring that a particle 
in a box by definition has ip = outside the box and that ip is everywhere continuous. See also section 3.2. 
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which may take any real non-negative value. There is also a lowest energy, given by setting 
n = 1. 



Definition When the possible energies of a quantum system are discrete and bounded below, 
the lowest possible energy is called the ground state energy (also sometimes called the zero point 
energy). The higher energies are, in increasing order, the first excited state energy, second excited 
state energy, etc. The corresponding wave functions are called the ground state wave function, 
kth excited state wave function. 

For the particle in a box the ground state energy is Ei = -k^Y? jlma^ , while E^ = Ei. Of 
course, we precisely wanted a theoretical understanding of such quantized energies in order to 
explain the energy levels of the hydrogen atom, determined empirically as (1.4). For the particle 
in a box, if we take m = m^- 2± 9.11 x 10~^^ kg to be the mass of an electron and a = 10~^° m 
to be the approximate size of an atom, we obtain 

En 0^ 5.97 X IQ-^^ra^ J . (2.22) 

In particular, the difference in energies between the ground state and first excited state is 
E2 — Eic^ 1.79 X 10~^^ J. A photon that is emitted in a transition between these energy levels 
then has a wavelength A ~ 1.12 x 10~^ m (on the boundary between the ultraviolet and X-ray 
parts of the electromagnetic spectrum), which is indeed the correct order of magnitude observed 
in atomic transitions! Here we are effectively modelling a hydrogen atom as an electron confined 
to an atom-sized box, which is very crude; we shall treat the hydrogen atom more precisely in 
section 5. Nevertheless, the above computation is encouraging. 
The full time-dependent wave functions (2.15) are 

r 5gi^mrxg-m27r2;it/2ma2 0<X<a, 

*n(x,i) = I - ' V, ■ (2.23) 

I U , otherwise . 

Any linear combination of such wave functions satisfies the time-dependent Schrodinger equation 
(2.13); in particular, the space of solutions, or possible wave functions, is infinite dimensional. 

In this example the energy levels E^, and associated stationary state wave functions '0^, are 
labelled naturally by a (positive) integer n. As we shall see throughout the course, such integers 
arise in many important solutions to the Schrodinger equation, and they are generally known 
as quantum numbers (although it is difficult to give a precise general definition). 

2.5 Degeneracy 

Having studied a particle in a one-dimensional box, it is now straightforward to extend this to 
a three-dimensional box. 

Consider a particle confined to the box region given hy {{x,y, z) e M.^ \ < x < a, < y < 
b, < z < c} C M.^ , where the potential is zero inside the box. In other words, 

0<x<a, 0<y<b, 0<z<c, 
V{x,y,z) = { . (2.24) 

otherwise . 
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As before, the stationary state wave function 'tp{r) = ijj{x, y, z) is taken to be zero on, and outside, 
the boundary of the box region. Inside the box the stationary state Schrodinger equation (2.12) 
reduces to 

+ Qy2 + 5^2 - ^2 ^ • y^-^^) 

This may be solved by separation of variables. Equation (2.25) essentially reduces to three copies 
of the one-dimensional equation, with the wave functions labelled by three quantum numbers 
ni, n2, na G 1^>o- Explicitly, inside the box these are given by 

Vm,n2,n3 {x, y,z) = B sm — — sm sm — — , (2.26) 

with B again an arbitrary constant, and the corresponding energies are 

Exercise (Problem Sheet 1) Derive the wave functions (2.26) and energies (2.27) by solving 
(2.25) by separation of variables. 

Definition If the space of solutions to the stationary state Schrodinger equation with energy 
E has dimension d > 1, wc say this energy level is d-fold degenerate; if it is one-dimensional we 
say E is a non- degenerate energy level. 

For the one-dimensional particle in a box all the energy levels are non-degenerate. However, 
consider now the three-dimensional box with equal length sides a = 6 = c, so that 

7^2 ^2 

In this case there are linearly independent wave functions with the same energy. For example, 
we may take (ni, n2, n^) to be any of (2, 1, 1), (1, 2, 1), (1, 1, 2), all of which have the same energy 
E = Gtt^ / 2ma'^ . There is thus a three-fold degeneracy in the number of quantum stationary 
states with this energy. The degeneracy in this case is related to the symmetry of the potential. 
We shall see other examples of this later. 

3 The Born interpretation 

We have now met the Schrodinger equation and solved it in the simplest interesting example, 
namely a particle confined to a box. We have seen that this leads to quantized energy levels, 
and that by crudely modelling a hydrogen atom as an electron confined to an atom-sized box, 
we obtain energies of the correct order of magnitude seen in atomic transitions. 

An immediate question is: what is the physical meaning of the wave function *(r,f) that we 
are solving for? Comparing to the corresponding classical problem, described at the end of 
section 2.1, we may also ask: where is the particle at time t7 In this section we shall answer 
these questions. 
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3.1 Probability density 

In order to motivate the interpretation that fohows, we begin by going back to the double slit 
experiment in section 1.3. In fact let us begin by discussing the corresponding experiment with 
light (also called Young's experiment). In this case the intensity of the light hitting the detector 
screen forms an interference pattern, and in classical electromagnetic theory this intensity is 
proportional to the absolute value squared of the amplitude of the wave. Unfortunately we won't 
have time to explain this in detail here, but this fact would have been well-known to the pioneers 
of quantum theory in the early 20th century. In the double slit experiment with electrons, we 
instead plot the spatial distribution of electrons hitting the detector screen over some long period 
of time, and then reinterpret this as a probability distribution for where any given electron will 
hit the screen. 

If we now conflate these observations, we are led to the hypothesis that the probability density 
function for an electron hitting the detector screen is given by the absolute value squared of the 
de Broglie wave associated to the electron. Of course, the detector screen could be anywhere, 
and the de Broglie wave is precisely the wave function ^(r,t) appearing in the Schrodinger 
equation. Hence we arrive at: 

Definition / postulate The function 



is a probability density function for the position of the particle, where ^'(r,t) is the particle's 
wave function. 

This interpretation of Schrodinger's wave function is due to Born. 

The assertion (3.1) is equivalent to the statement that the probability of finding the particle 
in a volume D C is given by 



This probability depends both on the region i), and also on the wave function ^'(r, t) satisfying 
the Schrodinger equation. We have already seen for the particle in a box that the space of 
solutions to the Schrodinger equation is an infinite-dimensional vector space, with basis (2.23). 
A given solution is said to describe the state of the particle, so that the probabilities (3.2) depend 
on the state. Notice P^{D) also depends in general on time t, although this is suppressed in the 
notation. 

Of course, the above assertions immediately raise some issues. In particular, the probability of 
finding the particle somewhere in should equal 1, at any time. Thus for (3.1), or equivalently 
(3.2), to make sense, the wave function must be normalized in the following sense: 

Definition A wave function ^ is said to be normalizable if 



p{v,t) = \^{v,t)\\ 



(3.1) 




(3.2) 




< oo , 



V time t . 



(3.3) 
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Moreover, if 

/ \^{r,t)fd^x = 1, Vtimei, (3.4) 
then ^' is said to be a normalized wave function. 

Similar definitions apply in one dimension. For example, the normalized condition (3.4) becomes 



/oo 
|*(x,i)pda; = 1. (3.5) 
-oo 



Let us make some remarks: 



1. The normalized condition (3.4) fixes the freedom to multiply a given solution to the 
Schrodinger equation by a complex constant, up to a constant phase e'*^. The latter 
is in fact not physical, so that wave functions differing by a constant phase are physically 
equivalent. Notice also that the normalized condition rules out the trivial solution = 0. 

2. For the particle in a box, the wave functions (2.23) are normalizable: 

r \^n{x,t)\''dx = \Bf Tsin^— dx = Ja|B|2. (3.6) 

The total wave functions 

, , , f . /f sin ^ e-in27r2ftt/2ma2 < x < a , , 

[ , otherwise , 

with B = y^2/a, are then normalized. Note that a particle being "confined to a box" may 
be interpreted as saying that there is zero probability of finding it anywhere outside the 
box, or in other words that the wave function is identically zero outside the box. 



3. More generally, notice that for a stationary state of energy E, we have 

p{r,t) = |^'(r,i)|2 = V'(r) e-^-^*/^ ^ = |V(r)P . (3.8) 



Thus *(r,t) is normalized for all time provided ip{r) is normalized, i.e. 

|V'(r)pd3x = 1. (3.9) 



/ 



/M3 

4. Recall that the plane wave 

n\k\H 



^{r,t) = A exp 



i k-r 



2m 



(3.10) 



satifies the free Schrodinger equation {V = 0), and was interpreted as describing a free 
particle of mass m and momentum p = - see the discussion around equation (2.3). 
However, clearly |^(r, t)\^ = |Ap is not integrable over M^, and thus the plane wave is not 
normalizable! In fact the plane wave is better interpreted as a beam of particles of mass 
m and momentum p. 
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* Some further, off-the-record, comments are in order about the plane wave. The plane 
wave is typically used in quantum mechanics in scattering problems. Here it is interpreted 
as an incoming beam of free particles of mass m and momentum p = ^k, which then scatter 
off some potential V that is non-constant only in a localized region of space. Readers who 
would like to learn more about this topic, which is elementary but beyond our syllabus, 
might like to read Chapter 5 of the book by Hannabuss. 

Alternatively, one can construct a localized (Gaussian) wave packet by taking an ap- 
propriate linear combination of plane waves. More precisely, this linear combination is 
realized as an integral over k. The resulting wave function satisfies the free Schrodinger 
equation and is normalized, although the wave packet spreads out as time increases. Again, 
we refer the interested reader to the references, rather than being more precise here. 

Finally, the real physical reason why the plane wave is not normalizable is because it 
describes a particle of definite momentum. In quantum theory one cannot simultaneously 
know the precise position and momentum of any particle, a fact known as Heisenberg's 
Uncertainty Principle; the more precisely one knows one quantity, the less precisely one 
knows the other. Since the plane wave particle has definite momentum, its position is 
completely uncertain, which explains why its formal probability density is a constant - 
it is equally likely be found anywhere in space. This fascinating implication of quantum 
mechanics is not something we will have time to explore in this course. This starred section 
is in any case getting ahead of ourselves, and the reader should perhaps ignore the above 
comments on a first reading, or indeed ignore them altogether if they are confusing. 



Example (Particle in a box) Let us apply the above ideas to the particle in a one-dimensional 
box. The normalized wave functions are (3.7). The corresponding probability density functions 
are then, for x G [0, a], given by 



Pnix) = - sin' 



mrx 1 2n7rx\ 

= - 1 — cos 

a \ o, J 



(3.11) 



a a 

and identically zero outside the box. Plots of the probability density functions for the ground 
state and first two excited states are shown in Figure 5. 

2 





Figure 5: Probability density functions for the first three states of a particle in a box. 

As usual in probability theory, we may define the distribution function as 
^ , , /"^ / N , X 1 . 2mrx 

Fn{x) = / pn{y)dy = - -—sin . 3.12 

Jq a zniT a 

Notice in particular that the first term, ^, is the result for the uniform distribution, where the 
particle is equally likely to be found anywhere in [0, a]. This is interpreted as the classical result. 
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By this we mean that if we are ignorant of the state of the classical particle before we observe 
its position, then because it moves at constant velocity back and forth across the box, it is 
equally likely to be found anywhere. The second, oscillatory term in (3.12) is then a quantum 
contribution, or quantum correction, to the classical result. 

Let us compute some example probabilities using these formulae. The probability that the 
particle is within | distance from the centre of the box is given by 



"|^„(x)pd. = fJ^-^)-F^[1 



1 1 , mr 

- + — sm — , 

2 mr 2 

, n even , 

1 (_l)(n-l)/2 

i + ' ^ odd . 



(3.13) 



In particular, we see that this approaches the classical result of |, for the uniform distribution, 
as n — > oo. The tendancy of quantum results to approach those of the corresponding classical 
theory for large quantum numbers is called the correspondence principle. 

We may similarly use the probability density p to compute expectation values: 

Definition In quantum mechanics the expectation value of a function of position /(r) is denoted 

(/(r))^ ^ / f{r)Mr,t)\'d^x. (3.14) 

The reason for the strange brackets on the left hand side will become apparent if you study 
quantum mechanics in the third year. Notice that again this expectation value depends on the 
wave function/state ^, and in general is a function of time t. 

Example For the particle in a box we compute the expected value of its position 

(^)*n = / XPnix)dx, 

Jo 

= [xFnix)]^- [ Fnix)dx , 
JO 



^x^ a 2n7rx"'" 
2a 4n^7r^ a 



= (3.15) 
agreeing with the classical result for the uniform distribution. 

3.2 The continuity equation and boundary conditions 

In this section we consider more carefully the boundary conditions involved in the Schrodinger 
equation. 

Proposition 3.1 The Schrodinger equation implies the continuity equation 

f^ + V-j = 0, (3.16) 
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where />(r,t) = |^'(r,f)p is the probability density we have already met, and the vector field 

i{r,t) = ^ (M/(r,i)(V*)(r,i) - *(r:t)(V^)(r,t)) , (3.17) 
is the probability current j. 
Proof This is a direct computation: 
^l*(r,t)|' = (|*M))^'(r,t) + *(M)|*(r,t), 



n \ 2m 



qi + qi 



i ( H 



2 



h V 2m 



, (Schrodinger) 



— V^^ + y* *-^- -— V^^' + F^ , (F is real) 
a \ 2m J n \ 2m ' 

]h 

2m ^ ' 
= ^V-(^V*-*V^) = -V-j. I (3.18) 

This leads to the following result: 

Proposition 3.2 Suppose thatM time t, j(r, t) satisfies the boundary condition that it tends to 
zero faster than l/|rp as |r| = r — )• oo. Then 

\^{r,t)\^d^x (3.19) 



/ 



is independent of t. In particular, if^ is normalized at some time t = t*, it is normalized Vt. 



Proof Let 5 be a closed surface that encloses a region D CM.^. Then 
d_ 
di 



f \'^{r,t)fd^x = f ^ d^x , (derivative through the integral) 

Jd Jd (yt 

= / (— V-j)d^x, (continuity equation) 
Jd 

j • dS , (Divergence Theorem) . (3.20) 



'S 

We refer the reader to appendix A for a reminder of the Divergence Theorem. In particular, 
dS = n dS" where d^ is the area element of S and n is the outward pointing unit normal vector 
to S. 

Now take 5 to be a sphere of radius r > 0, centred on the origin, so that D is a ball. Then 

n = - and hence 
r 

j-dS = j-nr^d^unit , (3.21) 

where d^unit is the area element on a unit radius sphere. In general, the function j • n will depend 
on r and the angular variables 6,(f) on the sphere (and on time t). Provided j • n = o(l/r^), 
uniformly in the angular coordinates, then using (3.21) the surface integral in (3.20) tends to 
zero as r — > oo, and hence 

1/ Mr,t)\'d^x = 0, (3.22) 
ot 7m3 

which implies that the expression in (3.19) is independent of t. | 
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We now state more formally the conditions that solutions to the Schrodinger equation should 
satisfy: 

1. The wave function ^{r,t) should be a continuous, single-valued function. This condition 
ensures that the probability density p = j^'p is single- valued and has no discontinuities. 
We already imposed the continuity property for the particle in a box. 

2. ^ should be normalizable, i.e. the integral of over all space should be finite. Propo- 
sition 3.2 ensures that if the probability current j in (3.17) tends to zero fast enough at 
infinity, then this integral is independent of time t. Then if ^ is normalizable, we may 
normalize it for all time. Without this condition, the probabilistic interpretation of the 
wave function is not possible. As already mentioned, this condition may be relaxed for 
free particles and scattering problems, although we shall not consider such wave functions 
any further in this course. 

3. should be continuous everywhere, except where there is an infinite discontinuity in the 
potential V . This condition follows since a finite discontinuity in V^' implies an infinite 
discontinuity in V'^^', and thus from the Schrodinger equation an infinite discontinuity in 
V. Again, we encountered precisely this behaviour for the particle in a box. 

3.3 Measurement of energy 

In section 2.2 we found separable solutions to the time-dependent Schrodinger equation (2.7). 
The corresponding stationary states (2.11) have definite frequency, and hence definite energy 
E. For example, the complete set of stationary state wave functions for a particle in a box is 
given by (3.7). Since the time-dependent Schrodinger equation is linear, any linear combination 
of such stationary state wave functions also solves the Schrodinger equation. Thus we may in 
general write 

*(r,t) = 5^a„Vn(r)e-^"*/\ (3.23) 

n 

where n labels some set of stationary states V'n(r) of energy En, and £ C are constants. 

For the particle in a box, (3.23) is simply a Fourier sine series. From Mods you know that 
any sufficiently well-behaved"^ function / : [0, a] — >• C with /(O) = /(a) = can be expanded as 
a Fourier sine series 

oo oo 

/(^) = ^ On a/ -sin = ^ani^nix) , (3.24) 

V 0- a 

n=l n=l 

for appropriate a„. Then given any such function f{x), we obtain a corresponding solution to 
the time-dependent Schrodinger equation using (3.23): 

oo 

^{x,t) = ^a„V'n(a;)e-'^"*/\ (3.25) 

n=l 

For the apphcation to wave functions below, we require the wave function to be noruializable as in (3.3). 
For the particle in a box example, this implies that f^ |/(a;)|^da: < oo. In fact this is then sufficient for the 
Fourier series (3.24) to converge almost everywhere to f{x). 
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where the coefficients an are determined via (3.24), and for the particle in a box the energies 
En are given by (2.21). We have thus solved the initial value problem for the Schrodinger 
equation, with '^{x,t = 0) = f{x). Note that the normalized stationary state wave functions 
fpnix) = sin are also orthonormal 

/ tpm{x) ll^n{x) dx = Smn ■ (3.26) 

Jo 

Definition / interpretation Suppose that the normalized wave function for a particle in a 
box is given by (3.25). Then the probability of measuring the energy of the particle to be En is 

This definition makes sense, since 

pa ^ pa ^ 

1 = / \^{x,t)fdx = V a^a„e-'(^"-^-)*/M ^„(x) Vn(x) dx = V |a„|2 . (3.27) 

•^0 „,,n=l Jo n=l 

Notice that if the particle's wave function is a stationary state of energy En, then the probability 
of measuring the particle's energy to be is 1. It thus makes sense to say that a stationary 
state (2.11) describes a particle of definite energy E. 

In fact the expansion (3.23) of any wave function '^{r,t) exists under very general conditions, 
i.e. for general choices of potential V{r), and generalizes to cases with degenerate energy 
levels or even a continuous energy spectrum. We have only discussed the quantum probabilistic 
interpretation of position and energy in this course, and our treatment of these two observables 
looks rather different. However, there is an elegant unified description of position and energy 
measurement that extends to all observables in quantum mechanics (including in particular 
momentum). This is described in the course B7.1a. 

4 The harmonic oscillator 

The quantum harmonic oscillator is ubiquitous in physics, and as such is probably the most 
important solution to the Schrodinger equation. 

4.1 The one-dimensional harmonic oscillator 

Consider a particle of mass m moving in one dimension under the influence of a potential V{x). 
Near to a critical point xq of V, where V'{xo) = 0, we have the Taylor expansion 

V{x) = V{xo) + lV"{xo)ix-xof + 0{{x-xof) . (4.1) 

If V"{xo) > then xo is a local minimum of the potential. Without loss of generality let us 
choose coordinates where the critical point xo is at the origin. Then to lowest order the potential 
near to xq = is 

V{x) = V{0) + ^mui^x'^ , (4.2) 
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where we have defined uP" = V"{0)/m. The dynamics is not affected by the value of V{0), since 
an additive constant drops out of the force F = —V'. The harmonic oscillator potential is 

V{x) = ^moj'^x^ . (4.3) 

The above analysis shows that any system near to a point of stable equilibrium is described by 
this potential. Classically, we have the force F = —V' = —mu^x so that Newton's equation 
(0.1) becomes 

x + J^x = , (4.4) 

which has general solution x{t) = Acosujt + Bsmoot. A classical particle in the potential (4.3) 
thus oscillates around the minimum with angular frequency uj. See Figure 6. 



VCx) 




^0 



Figure 6: A particle performing small oscillations around a point of stable equilibrium of any 
potential V{x) is described to lowest order by a harmonic oscillator. 



We now turn to the corresponding quantum mechanical problem. The stationary state 
Schrodinger equation (2.14) of energy E is 



d> 1 2 2, 



(4.5) 



The first thing to do is to redefine variables so as to remove the various physical constants: 

2E ^ I mco 



€ = 



h 



-X 



SO that (4.5) becomes 



where we have defined 



ex , 



moj 



h 



-X 



(4.6) 



(4.7) 



(4.8) 



Although at first sight (4.7) looks like a fairly simple ODE, depending on a single constant e, 
in fact it is not so simple to solve. However, it is not difficult to spot that xiO — e^^^^^ solve 
(4.7) with e = ±1. To see this, we compute 



d^ 



d^ 



(4.9) 
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As discussed in section 3, we are only interested in normalizahle solutions to the Schrodinger 
equation. For a stationary state, this means that 

/oo 
\ip{x)\^ dx < oo . (4.10) 
-oo 

Via the change of variable (4.6), (4.8), for x{C) = e^^^^^ the left hand side of (4.10) is 

jATeT^^dC, (4.11) 
V 7_oo 

which is finite only for the minus sign, giving a Gaussian integral. Thus the solution x{Cj = 
with e = —1, is not normalizahle, and we thus discard it. 
We now again change variables by defining 



X(0 = /(Oe- 



(4.12) 



so that 





fdl 


de 




d^X 


[dv 


de 





■d^ 



-e/2 



and hence (4.7) becomes 



dV 



24^ + (e-l)/ 







(4.13) 



(4.14) 



If anything, this looks worse than (4.7), so you might wonder why we bothered with (4.12)! 
The reasoning here is that for large values of |^| the term in (4.7) will dominate over the 
ex term; thus for large |^| one expects solutions for different e to have the same behaviour to 
leading order. This is why we have written xiO (4-12) as the above normalizahle e = +1 
solution times some other function /(C)- 

It is still not clear how to solve (4.14) exactly, so as usual for this type of equation we try a 
power series solution.^ Thus we write 



(4.15) 



fc=0 



and compute 



•dC 



fe=o 



and 



dV 

de 



(4.16) 



(4.17) 



fe=0 



* * This is a topic that could have been included in the Part A Differential Equations course, but wasn't. 

The point ^ = is an ordinary point of (4.14), and then Fuchs' theorem guarantees that any solution may be 
expressed as a convergent Taylor series (4.15). The interested reader is referred to appendix B. Having said this, 
most quantum mechanics textbooks suppress these details. 
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Notice that the last expression involves a relabelling of the original sum. By substituting (4.15), 
(4.16) and (4.17) into (4.14) we obtain 

oo 

^[(fc + l)(A; + 2)afc+2-2fcafe + (e-l)afc]e'= = 0. (4.18) 

A;=0 

The coefficient of every power of ^ must be separately zero, so we obtain the recurrence relation 

The even and odd powers are then decoupled, giving rise to the two linearly independent series 
solutions 



/even(0 = ^0 
foddiO = Ol 



■ , (l-€) ,2 , (5-€)(l-€) 4 

2! ^ 4! ^ 

f ^ .3 , (7-6)(3-e) ,5 

^ 3! ^ 5! ^ 



(4.20) 



Definition In one dimension a stationary state wave function satisfying ijj{—x) = ±'tp{x) is 

{even 1 
> parity state, 
odd J ^ 

Recalling equations (4.6), (4.8) and (4.12), the total stationary state wave functions are 
/even/odd(0 e^^^/^, which then describe even and odd parity states, respectively. 

We must now determine when these solutions are normalizable. It turns out that this is the 
case if and only if the series in (4.20) terminate. We shall not provide a full, rigorous proof 
of this here, but instead sketch the proof. This is purely for reasons of time. A more detailed 
treatment may be found in the book by Hannabuss. 

Notice first that either of the scries in (4.20) terminating is equivalent to the statement that 
/(^) is a polynomial, which in turn is equivalent to the statement that an+2 = for some integer 
n > 0, since then the recurrence relation (4.19) implies that = for all A; > n + 2. Note that 
ao and ai arc necessarily non-zero, otherwise /cvcn/odd('?) = 0) respectively. 

Suppose that either series in (4.20) does not terminate. Then all the coefficients are non-zero, 
and the ratio 

^ ^ ? , as ^ oo . (4.21) 

Ofc k 

Compare this asymptotic behaviour of the power series coefficients with that of the function e^^ : 



e^ = 



where we have defined the coefficients 



oo ^ oo 

EW^^'' = Jlhe , (4.22) 



q=0 ^' fe=0 



bk = I ^ " ' (4.23) 

lo, k = 2q+l. ^ ' 



For k = 2q even we then have 

bk+2 (I)! 2 2 



- yl^o^ - , ' , , as k ^ OO . (4.24) 

bk (^)! A; + 2 A;' ^ ^ 
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Comparing to (4.21), we thus see that /even(0 ^ h.ave the same asymptotic expansion. 
Hence the total stationary state wave function is 

xiO = /(C) e-«'/^ ~ e^"/\ (4.25) 

and thus ip{x) = xiO is not normaUzable.^ A similar argument applies to foddiO the 
function ^ e^^ . 

Thus for a normalizable solution the series in (4.20) must terminate. If n > is the least 
integer for which 0^+2 = then the recurrence relation (4.19) implies 

2n + 1 - e = . (4.26) 

Recalling the definition e = 2E/Huj in terms of the energy E in (4.6) then gives 

E = En = {n + l)nuj . (4.27) 

The quantum harmonic oscillator energies are hence labelled by the quantum number n G Z>o; 
compare to (2.21) for the particle in a box (and notice there that we instead defined n so that 
n G Z>o). Also note that our initial solution with / = 1 and e = +1 is in fact the ground state, 
with n = 0. Reverting back to the original spatial coordinate x via (4.6) and (4.8), the ground 
state wave function is hence 

V'o(x) = aoe-'"'^^'/^^ (4.28) 

where uq is a normalization constant, with corresponding ground state energy Eo = Ihu. The 
normalized ground state wave function may be obtained by imposing 

/oo 
Q-mu;xyh^^_ (4.29) 
-00 

The Gaussian integral is that for a normal distribution of variance a'^ = h/2mL0, and thus via 
the standard result for this integral we have 




(4.30) 



leading to the full, normalized time-dependent ground state wave function 

^o(x,i) = f"^y/'e-(--a.^+i/-t)/2/i. (4.31) 

More generally the stationary state wave functions are x„(C) = fn{0^~^^^^i where fn{0 
an even/odd polynomial in ^ = ^f^x of degree n, for n even/odd, respectively. The polyno- 
mials /„ may be determined explicitly by setting e — 1 = 2n in the recurrence relation (4.19). 
Appropriately normalized, /n(0 = -f^n(0 is called the nth Hermite polynomial. The first few 
polynomials are given in Table 1, with the corresponding wave functions shown in Figure 7. 

Exercise Check that the polynomials in Table 1 are indeed obtained from (4.20). 
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n 


1 

J. 


1 


2^ 


2 


4^2 _ 2 


3 





Table 1: The first four Hermite polynomials. 




Figure 7: The ground state and first three excited state wave functions of the harmonic oscillator. 



We now summarize via 



Theorem 4.1 The energies of the one- dimensional quantum harmonic oscillator of angular 
frequency oj are 

E = En= {n+\)huj , 

for n G Z>o a non-negative integer. The corresponding normalized stationary state wave func- 
tions are 

where ^ = \/mu/h x and is the nth Hermite polynomial, defined implicitly via the generating 
function: 

oo ^ 

exp [2^8 -s^) = HniO ^ • (4.33) 

n=0 

®It is the last step in (4.25) that needs a httle more work to make rigorous. 
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Exercise (Problem Sheet 2) You'll be led through the derivation of the normalized wave func- 
tions (4.32), and the expression (4.33) for the Hermite polynomials, on the second problem sheet. 
Note that for n = (4.32) agrees with the ground state result (4.31). 

Let us now discuss some physics. Notice that the difference in consecutive energy levels in 
(4.27) is En+i — En = hco, which is precisely the energy of a photon of angular frequency a;! 

* What is the relation between a quantum harmonic oscillator and a photon with the 
same angular frequency w? Recall that we said the photon is massless and therefore not 

governed by the Schrodinger equation. While this is certainly true, it turns out that 
in "quantizing" the classical electromagnetic field one effectively finds an infinite set of 
quantum harmonic oscillators. 

As stressed at the beginning of this section, the harmonic oscillator is ubiquitous in physics, 
and it is therefore straightforward to test these results experimentally. For example, in a diatomic 
molecule the two atoms perform small oscillations of some frequency io. One can then measure 
the frequency of absorbed/emitted photons, finding agreement with (4.27). 

As for the particle in a box, it is interesting to compare the classical and quantum harmonic 
oscillators. We focus here on the following new feature. In classical physics a particle of conserved 
energy E cannot enter a region of space where E < V{r), as follows from the simple observation 
that 

= ^ + V{r) > V{r) . (4.34) 

The set {r € | < ^(r)} is called the classically forbidden region, for fixed particle energy 
E. For example, the ground state of the quantum harmonic oscillator has energy Eq = ^Huj, so 
the corresponding classically forbidden region for this energy is {^huj < ^mu'^x'^}, which is the 
same as the set {|^| > 1}. The quantum probability of finding the particle in this region is 

/ \^o{x,t)\'^dx = ^/ e"^" d^ ~ 0.157 . (4.35) 

mu) -'ICI>1 

Here wc have used (4.31) the the harmoinic oscillator ground state wave function. We thus 
sec that there is a non-zero probability of finding the particle in a region where classically it is 
impossible to find the particle! This is what underlies the phenomenon of quantum tunnelling, 
although wc shall not pursue this further in this course (regrettably). In keeping with the 
correspondence principle, for the analogous computation with energy En one finds the probability 
of finding the particle in the classically forbidden region tends to zero as n oo. We shall not 
prove this though. 

4.2 Higher dimensional oscillators 

Having discussed the one-dimensional oscillator, it is now straightforward to extend this to 
oscillators in higher dimensions. 

Let us begin by considering a quantum harmonic oscillator in two dimensions, with potential 

V{x,y) = ^m{oJix''+ojly'') . (4.36) 
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The corresponding stationary state Schrodinger equation (2.12) is 



(^'t + »^]+V(.,y)^ = E^. (4.37) 



2m \ dx^ dip . 

As for the particle in a box, this may be solved by separation of variables. One writes il){x, y) = 
X{x)Y{y), so that (4.37) separates into 



(PX 1 



+ -mulx^X = EiX , 



2m, dx^ 2 



where E\ ^ = E. We thus have two decoupled one-dimensional oscillators, and it follows 
from the previous section that 

E = En^^ri2 = [ni + l)nLOi + (n2 + I)nuj2 , (4.39) 

where the quantum numbers ni, n2 G Z>o- The corresponding normalized stationary state wave 
functions are 



ikm,n2{x,y) = ipni{x) Ai2iy) (4-40) 

where ipn denotes the normalized stationary state wave function for a one-dimensional harmonic 
oscillator. 

Exercise Check carefully that the above statements follow from Theorem 4.1. In particular, 
how is the normalizable/normalized condition related to the condition in one dimension? 

Of course the quadratic form appearing in (4.36) is rather special, in that it is diagonal. One 
can treat more general quadratic forms by first changing to normal coordinates, i.e. one first 
diagonalizes V by an orthogonal transformation. Let us illustrate with an example: 

Example Consider the oscillator potential 

V{x,y) = muj^ix'^ +xy + y'^) = ^mw'^i^x y)i^ \ \ ^ [j^y^ ' ^^'^^^ 

The matrix here has eigenvalues 3 and 1, so there exists an orthogonal transformation to coor- 
dinates u, V with corresponding new potential 



u V 



3 \ u 
I I \ V 



= {Su'^u^ + oj^v"^) , (4.42) 

which is now of the form (4.36).^ The orthogonality of the transformation ensures that the 
Laplacian/kinetic term in the Schrodinger equation is invariant, so 



^Explicitly, u ^ ^{x + y), v = ^{x - y). 
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Of course, since any quadratic form for V is described by a symmetric matrix, we may always 
diagonalize the potential by an orthogonal transformation, thus reducing the problem to a 
diagonal form for V, as in (4.36). These remarks apply in any dimension: first change to normal 
coordinates, then separate variables one at a time. The Schrodinger equation for any quadratic 
potential then reduces to a decoupled set of one-dimensional oscillators, so that the total energy 
= sum of one-dimensional energies, and the total wave functions = product of one-dimensional 
wave functions. 

We conclude this section with another example of degenerate energy levels: 

Example (Degeneracy) Consider the original two-dimensional oscillator potential (4.36) with 
ui = UJ2 = oj. Then 

E = En = in + l)nLU , (4.44) 

where n = ni + n2- The ground state is ni = 77,2 = 0, which is the unique state of energy hw. 
However, more generally at level n there are n + 1 linearly independent wave functions with 
energy En, given by taking (ni,n2) to be (n, 0), (n — 1, 1), . . . , (l,n — 1), (0,n). Thus En has 
degeneracy n + 1. Again, notice this degeneracy is related to the symmetry of the potential. 

5 The hydrogen atom 

5.1 Atoms 

We begin with a discussion of Coulomb's law J In general, a point charge 62 at the origin induces 
an electrostatic force on another point charge ei at position r given by the inverse square law 

F = 2" - ' 5.1 

where as usual r = |r|. Notice that the Coulomb force (5.1) is proportional to the product of the 
charges, so that opposite (different sign) charges attract, while like (same sign) charges repel. 
Electric charge is measured in Coulombs, C, and the proportionality constant eo — 8.85 x 10~^^ 
C^ m~^ in (5.1) is called the permittivity of free space. The Coulomb force is conservative, 
of the form F = — where 

= = 4^ • <^-^' 

is the Coulomb potential. 

An atom consists of negatively charged electrons e~ orbiting a positively charged nucleus, 
where the force of attraction is (predominantly) electrostatic. An electron carries electric charge 
— e, where e ~ 1.60 x 10^^^ C. In general, the nucleus of an atom consists of Z positively 
charged protons, each of charge +e, and A neutrons, of charge 0, and these are tightly bound 
together (by the strong force). In the following we consider an atom of atomic number Z, with 

* Coulomb discovered his law in 1783. From a modern point of view, this law is subsumed into Maxwell's 
classical theory of electromagnetism (i.e. Coulomb's law can be derived from Maxwell's equations). 
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a single orbiting electron. In particular, the hydrogen atom has Z = 1. Since the nucleus 
is vastly heavier than the electron (mproton — 1836 mg-), this implies that the centre of mass 
of the atomic system will always be very close to the nucleus. We thus make the simplifying 
assumption that the nucleus is fixed, at the origin. The potential for the force acting on the 
electron is then given by (5.2) with ei = — e, 62 = Ze (see Figure 8): 

Ze^ 
47reor 

sing/e electron 

^ charge = -e 



, Coulomb force 



lucleus: ' 

Z protons 
\A neutrons , 

charge 62 = Ze 

Figure 8: An atom consisting of a nucleus of atomic number Z and a single orbiting electron. 

Recall that in section 1.2 we discussed the empirical formula (1.4) for the energy levels of the 
hydrogen atom. Our task in the remainder of this course is to derive this formula theoretically, 
using the Schrodinger equation. However, before moving on to this, as usual we pause to make 
some comments on the corresponding classical problem. The 1/r Coulomb potential is formally 
the same as Newton's gravitational potential that you encountered in Mods, the only difference 
being that gravity is always attractive. Thus classically the atomic system in Figure 8 is similar 
to a planet orbiting the sun under gravity.^ In particular, the energy levels are continuous, not 
discrete as in (1.4). 

* There is an even more serious problem with the classical picture above. Maxwell's clas- 
sical theory of electromagnetism predicts that an accelerating charge, such as an electron 
in orbit around a nucleus, emits electromagnetic radiation, and thus continuously loses 
energy. The electron would quickly spiral in towards the nucleus, and all atoms would 
hence be unstable! We shall see how quantum theory gets around this problem. 

5.2 The spectrum of the hydrogen atom 

We consider the three-dimensional stationary state Schrodinger equation (2.12) with potential 
(5.3): 

-;^VV-/^V' = E^, (5.4) 
2m 47reo r 



Of course there is also a gravitational attraction between the electron and nucleus, but this is many orders 
of magnitude smaller. 
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where m = m^- ~ 9.12 x 10~^^ kg is the mass of an electron. In this course we consider only 
spherically symmetric wave functions, so that '0(r) = il){r) depends only on the distance r from 
the origin.^ In this case, the Laplacian reduces to 

W = -—iri.) = ^ + 2dV^, (5.5) 
r dr^ dr^ r dr ' 

and we must hence solve the ODE 

1 d2 ^ 2mZe2 , 2mE , , . 

rd^^^^^^ + iW^^ = ^^-^^ 

Notice that by multiplying this equation by r we effectively obtain a one- dimensional Schrddinger 
equation of the form (2.14), with ^p replaced by rtp and a 1/r potential. 

The method for solving (5.6) is very similar to that for the harmonic oscillator - there is 
essentially only one technical difference, as we shall see. As for the harmonic oscillator, we first 
change variables to remove the physical constants. We define 

where we have introduced 
Definition The Bohr radius is 



2 



a = ^ 

where m = m„- is the mass of the electron. 



2± 5.29 X 10-^ V , (5.8) 



Of course, a priori the constant k may be imaginary, but we are anticipating that E will be 
negative, as in (1.4). Equation (5.6) now becomes 

-r^ + -^ + -'^ = K tp . (5.9) 
dr"' r dr r 

As for the harmonic oscillator, it turns out to be convenient to substitute 

V'(r) = /(r)e— . (5.10) 

The motivation is again similar.^'' At large values of r the term on the right hand side of 
(5.9) dominates over the second and third terms on the left hand side. One thus expects the 
■0(r) = A±e^'^^ solutions to the equation d'^ip/dr'^ = K^ip to determine the leading asymptotic 
behaviour of solutions to (5.9). The normalization condition is 

r poo 

1 = / \'4}{r)\^d^x = \ip{r)\'^ r^dr , (5.11) 

JRS Jo 

where we have taken ij;{r) = ipi^), a'Hd integrated over the unit sphere of area 47r. Thus one 

must cortaiuh' take thc^ minus si"n in e^'^'' to obtain a normalizable solution near to r = oo. 



^The more general case is studied in the third year course 
■'The argument in the following paragraph is again re 
any case simply a justification for the substitution (5.10) 



^°The argument in the following paragraph is again really a sketch, rather than a rigorous treament. It is in 
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We also see that E must be negative; otherwise |e^'*''| = 1, and this asymptotic solution is not 
nor mahz able at infinity. 

The substitution (5.10) leads to 

-Kf]e 



so that (5.9) becomes 



dip 




dr 


\dr 


dV 




dr^ 


\dr^ 



+ ) e-'^'' , (5.12) 



d^f 2 df 1 

^ + -(l-.r)^ + -(/?- 2.)/ = 0. (5.13) 

We again try a series solution. The main technical difference with the harmonic oscillator 
equation is that the coefficients in (5.13) are singular at r = 0.^^ Thus a priori we should write 

oo 

fir) = i;ofcr^+% (5.14) 
fc=o 

for some real constant c G M. Without loss of generality, we may assume oq 7^ 0. We then have 

df 



dr 
dr^ 



5^(fc + c)afer^+^-i , 



A:=0 
2- °° 



^ -- Y.ik + c){k + c-l)akr''+^-'' . (5.15) 



k=0 

Thus (5.13) becomes 



{ik + c){k + c — 1) Ofe + 2[{k + c) Ofe — K{k + c — 1) 0^-1] 

fc=0 

+(/3-2/«)afe_i}r*^+^-2 = , (5.16) 
where we define a_i = 0. The lowest power of r in (5.16) is r'^~^, with coefficient 

ao[c(c-l) + 2c] = 0. (5.17) 

Since oq / 0, we thus have 

c(c-l) + 2c = ^ c(c+l) = ^ c = or c = -1 . (5.18) 

However, if c = — 1 then from (5.14) the wave function tfj{r) = f{r) ~ ^ near to r = 0. The 
wave function is singular at the origin with this choice of c, and we thus discard this solution. 
Hence c = 0, and the recurrence relation derived from the coefficient of r^~'^ in (5.16) simplifies 
to 

k{k-l)ak + 2[kak- K{k-l)ak-i] + {P -2K)ak-i = 0, (5.19) 



* More precisely, r = is a regular singular point of (5.13). Again, Fuchs' theorem applies and determines 
the form of the series in (5.14). See appendix B. 
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or, after some cancellations, equivalently 



2Kk- P , , 



The solution is hence 



f(.r) = ao 



' (2k -/3) (4k-/3)(2k-/3) 2 
1 + ^ — ^ r + '-^ — + 



(5.21) 



Notice that, in contrast to the harmonic oscillator, here we have only a single solution - we have 
already discarded the c = — 1 solution as singular and thus unphysical. 

Finally, we must impose that the wave function resulting from (5.21) is normalizable. Suppose 
that the series does not terminate. Then 

— > , as A; ^ 00 . (5.22) 

Compare this to the Taylor expansion of the function e^**^ = X^feLo^fc^'^' where 5^ = {2K)^/k\. 
In this case 

h _ {2Kf{k-l)\ _ 2_K 

bk-1 ~ (2K)fc-iit! " k ■ '^^ 

Thus if the series (5.21) does not terminate the function /(r) and e'^'^^ have the same asymptotic 
expansion, and hence 

V,(r) = /(r)e-'^^ ~ e"'' (5.24) 

is not normalizable in (5.11). 

The series must hence terminate; that is, there is a least integer n > such that a„ = 0. The 
recurrence relation (5.20) then gives 

2Kn-P = . (5.25) 
Recalling the definitions (5.7) of k and (3, this becomes 

so that 

" 2ma^ ' ^ ^ 

which is precisely of the form (1.4)! Moreover, we have a theoretical formula for the Rydberg 
constant, namely (putting Z = 1 for hydrogen) 

^ " 47r mca2 " {Anfelh^c ' ^^"^^^ 
where m is the electron mass, c is the speed of light, and a is the Bohr radius. 

In general /(r) = /n(r) is a polynomial of degree n — 1. In particular, the ground state has 
n = 1 with /(r) = oq a constant. The corresponding ground state wave function is 

^^{r) = aoe-'^i^ = 00 6-^"^/" , (5.29) 
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where we have used (5.25) to write k = Kn = P/2n = Z/an. It is straightforward to normahze 
(5.29) using (5.11): 



47r |ao| 



r^e-^'^i'^dr , 







d2 
2 CI 


TT 


ao 








2 


TT 


ao 


dKl 








TT 







-^''i'- dr , 





1 

2ki 



Z3 



thus giving oq = ^JT/^Jtuj? and hence 

■01 (r) = 

is the normalized ground state wave function. 



-Zr/a 



(5.30) 



(5.31) 




Figure 9: The ground state and first three excited state wave functions of the (spherically 
symmetric) hydrogen atom. 

More generally, the polynomials fn, appropriately normalized, are called generalized Laguerre 
polynomials. More specifically, 



(2Zr\ 
—) = Li_^ (S/^nr) , 

where the first few polynomials L\_-^{s) are shown in Table 2. 

Exercise Check that the polynomials in Table 2 are indeed obtained from (5.21) 



(5.32) 
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n 




1 


1 

± 


2 


-s + 2 


3 


^ - 3s + 3 


4 


+ 2s2-6s + 4 



Table 2: The generalized Laguerre polynomials that determine the ground state and first three 
excited states of the (spherically symmetric) hydrogen atom. 



* The generalized Laguerre polynomials may be defined via 

This is a polynomial of degree N. Our case above is = n — 1, a = 1, although we shall 
not prove that (5.33) follows from (5.21) here. Related to this, you might wonder how 
general the spherically symmetric case is. Allowing = ip{r, 9, 0) in (5.4) one finds that 
the spectrum of energies (5.27) is the same, but now En has degeneracy v?. The general 
wave functions involve spherical harmonics and the more general Laguerre polynomials in 
(5.33), with a > 1. The proof of these facts is in the third year quantum course. 

We summarize with 

Theorem 5.1 The (spherically symmetric) energy levels for an atom consisting of a single 
electron orbiting a nucleus of atomic number Z are 

where the positive integer n is called the principal quantum number. The corresponding spheri- 
cally symmetric wave functions are 

Mr) = C„Li_i(^)e-^'-/«", (5.35) 

where Cn are appropriate normalization constants and L\_-^ is a polynomial of degree n — 1. 
In particular, the normalized ground state wave function is 



Z3 



Mr) = \/^e-^^/V (5.36) 



Notice that there is a minimum energy ground state, with tt, = 1, which leads to a stable 
atom. We conclude with an example computation using the hydrogen atom ground state wave 
function (5.36): 

Example Let us compute the expectation value of r, the distance of the electron from the 
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nucleus, in the ground state. This is given by 

/•oo 
^0 



_47rag T 1 
8^ d^f |.2kiJ ' 

'24 ' 

where we have used oq = \/^^V^ra^ in (5.29). Thus for the hydrogen atom, with Z = 1, the 
average distance of the electron from the nucleus, in the ground state, is | times the Bohr radius 
a, which numerically is 2± 7.94 x lO"^-*^ m. 

A Summary: vectors and vector operators 

The following summarizes some basic notation and conventions used in the main text. 

We work in M^, or a domain therein, in Cartesian coordinates. If ei = i = (1, 0, 0), 62 = j = 
(0, 1, 0), 63 = k = (0, 0, 1) denote the standard orthonormal basis vectors, then a position vector 
is 

3 

r = ^XiCj , (A.l) 

i=l 

where xi = x, X2 = y, X3 = z are the Cartesian coordinates in this basis. The scalar product of 
two vectors a, b is denoted by 

3 

a ■ b = ^ aibi . (A.2) 

i=l 

We denote the Euclidean length of r by 

|r| = r = ^Jxj + xl + xl = Vr • r . (A.3) 
The gradient of a function (p = 0(r) is the vector field 

grade/) = V0 = ^^e^. (A.4) 

i=l ' 

The divergence of a vector field f = f(r) is the function (scalar field) 

divf = V.f = i:e « =i:|^. (A.5) 



dxi ^ dxi 

i=\ i=l 



The second order operator defined by 



W = V.(V^) = (A.6) 

i=l * 



is called the Laplacian. 

We also remind the reader of 
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Theorem A.l (Divergence Theorem) Let D be a bounded region in with boundary surface 
S = dR. If f is a vector field then 

' V-fd^x = / f -nd^ (A.7) 

ID Js 

where n is the outward unit normal vector to S = dR. One sometimes writes dS = ndS. 

B * Fuchs' theorem and Frobenius series 

This appendix is included only for interest and completeness. It is certainly not intended to be 
part of the course. 

In the main text we solved both the harmonic oscillator (4.14) and spherically symmetric 
hydrogen atom (5.13) ODEs using a (generalized) power series expansion. Here we discuss this 
problem more generally. Consider the second order ODE 

= 0, (B.1) 

for the function f{x). A point is called an ordinary point of (B.l) if p(a;), q{x) have Taylor 
expansions about xq, valid in some positive radius of convergence; otherwise xq is called a 
singular point. If xq is a singular point, but (x — xq)p{x), (x — xo)^g(x) both have Taylor 
expansions around xq, then xq is called a regular singular point. 

By a generalized power series, or Frobenius series, about Xq we mean 

oo 

/(x) = 5^afc(x-xo)'=+% (B.2) 

fc=0 

for some real number c G M. Without loss of generality, ao 7^ 0. Notice that for c a non-negative 
integer this reduces to a normal power series. Without loss of generality, we now set xq = 0. The 
equation obtained from the ODE (B.l) by setting the lowest power of x to zero in the expansion 
(B.2) is called the indicial equation. Since (B.l) is second order, this is a quadratic equation for 
c. Fuchs and Frobenius proved the following result: 

Theorem B.l 

• If X = is an ordinary point, then (B.l) possesses two distinct power series solutions, so 
that one can effectively set c = in (B.2). These series converge for \x\ < r, where r is 
the minimum radius of convergence of p{x) and q{x). 

• If x = is a regular singular point, then there is at least one solution of (B.l) of the form 
(B.2). Again, this series solution converges for |x| < r, where r is the minimum radius of 
convergence of xp{x) and x^q{x). 

The precise behaviour for regular singular points is a little involved. If the indicial equation 
for c has distinct roots that differ by a non-integer, then there are two series solutions of the 
form (B.2), with the corresponding values of c solving the indicial equation. Otherwise there is 
not necessarily a second series solution; if there is not one can nevertheless say more about the 
second solution, but we refer the interested reader to the literature for details. 
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